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摘 要 : 基于 1971 一 2015 年 陕西 省 无 缺 测 情况 分布 较 均匀 的 58 个 气象 站 点 的 逐日 降水 数据 ,采用 最 大 1 日 .3 
日 5 日 和 7 日 的 降水 量 代表 极端 降水 ,应 用 地 区 线性 矩 法 研究 区 域 极 端 降水 的 时 空 特 征 。 结 果 表 明 :(1) 陕西 省 可 
划分 为 6 个 水 文 气象 一 致 区 ,其 中 广义 极 值 分 布 (GEV ) 在 各 一 致 区 的 模拟 效果 最 好 ,各 一 致 区 最 优 频率 估计 值 与 同 
频率 的 实测 值 较为 吻合 。(2) 地 区 分 析 法 计算 的 极端 降水 频率 估计 值 较 单 站 分 析 法 具有 更 好 的 稳健 性 和 准确 性 ,万 
其 是 在 计算 较 长 时 段 的 极端 降水 情况 下 更 显著 。(3) 重 现 期 为 2a 一 遇 时 , 陕 南 的 地 区 增长 因子 大 于 陕 北 ; 重 现 期 达 
到 5 a 一 遇 时 则 相反 , 且 随 着 重 现 期 的 增长 ,地 区 增长 因子 及 其 在 陕 南 、 陕 北 的 差异 也 在 增 大 ,(4) 在 100a 和 50 a 一 
遇 重 现 期 下 ,极端 降水 在 陕西 省 南部 较 大 ,东部 居中 ,中 部 的 咸阳 至 商洛 地 区 西部 的 延安 西北 和 榆林 西部 较 小 ; 极 
端 降水 的 分 布 特征 与 陕西 省 独特 的 地 理 特征 有 关 , 尤 其 是 东西 向 分 布 的 秦岭 阻挡 了 水 汽 向 北 输送 ,造成 了 极端 降 


水 的 南北 差异 。 
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在 全 球 气候 变 暖 背景 下 ,极端 降水 事件 发 生 的 
频次 明显 增加 ,对 农业 生产 .社会 经 济 .生命 安全 
和 生态 环境 等 方面 造成 了 严重 的 影响 。 近 年 来 , 针 
对 区 域 极端 降水 的 现状 , 翟 盘 成 等 “通过 对 中 国 近 
几 十 年 极端 降水 事件 特征 的 分 析 ,指出 极端 降水 的 
平均 强度 和 极端 降水 事件 的 发 生 频 率 都 呈 增 大 趋 
势 。 对 于 地 区 性 极端 降水 的 特征 ,众多 学 者 探索 并 
分 析 其 原因 ,如 杨 涵 消 等 "研究 发 现 , 西 太平 洋 暖 池 
异常 升温 造成 的 海陆 温差 减 小 ,是 中 国 东 部 夏季 极 
端 降水 1990 年 前 后 出 现 转变 的 重要 驱动 因素 之 
一 。 龙 妍 妍 等 ”指出 ,高原 夏 季风 的 强 弱 会 造成 我 
国 不 同 地 区 夏季 极端 降水 的 增 减 , 而 中 高 层 环流 系 
统 特征 、 中 低层 风 场 特征 、 水 汽 输送 特征 和 湿 位 涡 
特征 是 其 对 极端 降水 影响 的 物理 机 制 。 

亦 有 大 量 学 者 对 陕西 省 极端 降水 事件 开展 研 
究 , 如 郑 小 华 等 "发 现 ,陕西 年 最 大 降水 量 和 强 降水 
发 生 频次 均 呈 现 增 加 趋势 ,并 分 别 在 20 世 纪 70 年 
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代 中 期 和 80 年 代 中 期 发 生 突 变 ; 王 雯 燕 '" 发 现 , 陕 
西 汛期 极端 降水 事件 发 生 频次 的 第 一 模 态 是 北部 
与 其 他 地 区 呈 反 位 相 变化 ,第 二 模 态 是 以 秦岭 为 界 
的 南北 间 反 向 变化 。 李 双双 等 "分析 了 秦岭 地 区 4 
种 极端 降水 类 型 的 时 空 变化 特征 ,进而 指出 当 气 温 
升 高 时 ,持续 性 极端 降水 呈 下 降 趋势 , 单 日 型 极端 
降水 呈 增 加 趋势 。 肖 科 丽 等 所 则 从 环流 与 海 温 蜡 
常 分 析 了 2010 年 陕西 夏季 极端 多 十 的 原因 。 可 见 ， 
陕西 省 极端 降水 的 研究 中 ,气象 学 者 从 极端 降水 的 
类 型 .强度 .频次 和 成 因 等 方面 进行 了 深入 人 研究 ,得 
出 了 一 些 重要 的 结论 。 

陕西 省 地 处 我 国内 陆 腹 部 ,南部 以 山地 为 主 ， 
北部 以 黄土 高 原 为 主 , 横 跨 湿润 - 半 湿 润 区 和 干旱 - 
半 干 旱 区 ,其 气候 特征 具有 过 渡 性 ,是 研究 气候 变 
化 的 典型 区 域 。 研 究 陕西 省 的 极端 降水 特征 ,可 以 
为 全 球 变 暖 及 水 循环 变异 特征 提供 理论 基础 。 针 
对 极端 降水 的 频率 分 析 ,通常 是 对 其 进行 线性 拟 合 
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和 参数 估计 等 ,进而 得 出 极端 降水 频率 估计 值 , 其 
极端 降水 估计 方法 包括 常规 矩 、 概 率 权重 矩 线性 
和 矩 和 地 区 分 析 法 等 。 地 区 线性 矩 法 是 在 线性 矩 的 
基础 上 ,对 研究 区 进行 了 一 致 区 划分 ,并 进一步 佑 
计 站 点 的 降水 频率 设计 值 ,是 当前 广泛 应 用 的 区 域 
极端 降水 频率 分 析 方 法 ,可 显著 提升 参数 估计 的 精 
准 性 和 频率 佑 计 的 准确 性 ,同时 ,水 文 气象 一 致 
区 的 划分 中 综合 考虑 了 水 文 气象 和 地 理 条 件 (经 纬 
度 、 高 程 和 地 形 地 貌 等 )。 该 方法 已 被 国内 不 少 学 
者 用 于 地 区 极端 降水 的 分 析 ”"”。 目 前 ,针对 陕西 
省 进行 极端 降水 的 频率 分 析 中 , 牛 晨 ”基于 线性 矩 
法 ,选用 3 种 线 型 对 关中 平原 极端 降水 的 频率 进行 
了 深入 分 析 ,其 频率 佑 计 值 与 实测 值 有 较 好 的 一 臻 
性 。 但 对 于 整个 陕西 省 ,还 没有 相应 的 极端 降水 地 
区 频率 人 研究。 因此 ,本 研究 基于 研究 区 新 近 的 逐日 
降水 资料 ,采用 地 区 线性 矩 法 分 析 陕 西 省 极端 降水 
的 重 现 期 和 时 空 分 布 特征 ,可 综合 分 析 陕 西 省 不 同 
气候 区 的 极端 降水 频率 ,对 于 区 域 防 灾 减 灾 具 有 一 


1 数据 与 方法 


1.1 数据 

本 文选 取 陕 西 省 58 个 气象 站 点 的 降水 数据 与 
陕西 省 的 DEM 资料 ,气象 数据 为 国家 气象 信息 中 心 
提供 的 逐日 降水 资料 ,序列 长 度 为 1971 一 2015 年 , 气 
象 站 点 分 布 见 图 1。DEM 资 料 来 自 地 理 空间 数据 云 。 

利用 陕西 省 1971—2015 年 共 45 a 的 逐日 降水 
资料 ,提取 出 58 个 站 点 的 年 最 大 1 日 降水 (AMP1)、 
年 最 大 3 日 降水 (AMP3) ,年 最 大 5 日 降水 (AMP5) 
和 年 最 大 7 日 降水 (AMP7) ,并 基于 此 研究 陕西 省 极 
端 降水 时 空 特征 。 
1.2 研究 方法 
1.2.1 地 区 线性 手法， 地 区 线性 矩 法 包括 线性 矩 法 
和 地 区 分 析 法 ,其 中 ,线性 矩 法 是 一 种 新 的 参数 佑 
计 方 法 ,其 估计 参数 具有 很 好 的 无 偏 性 和 稳健 性 ; 
而 地 ee 区 内 所 有 站 点 的 历史 降水 
资料 ,其 基本 思想 是 “以 空间 来 换 时 间 ” ,可 以 克服 
E ,通过 分 析 一 致 区 内 
气象 站 点 的 频率 分 布 曲 线 ,进而 计算 出 各 站 点 的 降 
水 频率 估计 值 。 地 区 线性 矩 法 综合 了 上 述 2 种 方法 
的 优点 ,已 经 在 国内 外 很 多 地 区 成 功 运 用 ”” i。 该 
方法 的 详细 介绍 ,参考 Hosking 等 的 研究 |。 


海拔 高 度 /m 
3748 
| 125 


a 气象 站 点 


>z 


图 1 陕西 省 58 个 气象 站 点 的 空间 分 布 示 意图 


Fig.1 Spatial distribution of 58 meteorological stations 


in Shaanxi Province 


线性 矩 的 概念 由 Hosking2 提出 ,其 基本 思想 是 
计算 服从 某 种 分 布 函数 的 次 序 统计 量 的 期 望 值 ,并 
对 其 进行 线性 组 合 。 假 设 世 服从 某 一 分 布 函数 ， 
Xin Xon < Xp _ 是 次 序 统计 量 , 将 > 阶 线性 矩 变 量 
定义 为 : 

A, =r S ('C EX, u r=1,2, n (1) 


式 中 : X e 是 样本 容量 为 x 、. 排 在 第 -位 的 次 序 

统计 量 ; EX, py 为 次 序 统 计量 X, p 的 期 望 值 ; Ch, 

JEMA r- 1 个 不 相同 的 元 素 中 选取 大 个 元 素 的 组 合 数 。 
区 的 前 四 阶 线性 矩 为 : 


fees i 
li=n 


=j+l j=l (2) 
= 2X a Xia) 


-1 2 3 


l= 上 5 a = 


式 中 : 为 样本 容量 ,基于 依次 增加 的 次 序 

X <X <…, ,为 样本 次 序 统 计量 ;利用 前 四 阶 
样本 线性 矩 可 推 求 其 统计 特征 参数 ,包含 离 差 系数 
(L-C, ) 、 偏 态 系数 ( 工 -C,) 和 峰 度 系数 ( 工 - C，)， 
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ae ee 1.2.3 最 优 分 布线 型 的 确定 ”水文 领域 上 通常 使 用 

l 2~5 个 参数 的 分 布线 型 进行 频率 计算 ,其 中 ,三 参数 

a (3) 分 布线 型 在 模拟 极端 降水 的 分 布 时 不 仅 准确 ,而 且 

i 应 用 灵活 方便 。 在 地 区 线性 矩 方法 中 , 三 参数 分 布 

L-C,=t=7 线 型 中 的 广义 逻辑 斯 蒂 分 布 (GLO) ,广义 极 值 分 布 


地 区 分 析 法 的 研究 对 象 必 须 属于 水 文 气象 一 
致 区 ,其 在 极端 降水 特性 上 具有 相同 的 气候 背景 和 
统计 特性 。 该 方法 将 降水 序列 分 解 为 两 部 分 , 即 反 
映 本 站 点 降水 特性 的 本 地 分 量 和 反映 区 域 降水 共 
性 的 区 域 分 量 。 推 算 站 点 的 频率 估计 值 时 ,综合 
虑 这 2 种 分 量 , 使 得 计算 结果 更 加 可 信 。 频 率 估 计 
值 的 计算 公式 如 下 : 

Qnij= 0X; (4) 
式 中 :7 为 所 要 推 求 极端 降水 的 重 现 期 ;i 代表 各 一 
致 区 ; j 为 相应 一 致 区 内 的 站 点 ; gy 为 一 致 区 ;的 
频率 因子 ,反映 该 区 各 站 点 共有 的 极端 降水 特性 ， 
称 为 地 区 频率 因子 或 地 区 增长 因子 ,可 利用 该 区 的 
最 优 分 布 曲线 和 分 布 参数 求 得 ,而 分 布 参数 可 通过 
区 域 线性 矩 系数 推 求 ; x 为 本 地 分 量 因 子 ,为 i 区 j 
站 点 极端 降水 量 的 多 年 平均 值 。 根 据 公 式 (4) 
可 得 : 


pa (5) 
即 站 点 频率 估计 值 经 过 去 均值 化 ,就 可 以 得 到 
反映 地 区 共性 的 极端 降水 分 量 , 即 地 区 增长 因子 。 
1.2.2 水 文 气象 一 致 区 的 划分 ”为 使 一 致 区 内 的 站 
点 具有 相同 的 水 文 气象 背景 ,本 文 从 以 下 3 个 方面 
对 站 点 进行 一 致 区 的 划分 : 
(1) 气象 学 上 的 判定 , 即 一 致 区 内 的 各 站 点 在 
水 汽 和 人流 和 物理 成 因 等 背景 条 件 上 要 相同 ,这 方面 
通常 借助 海拔 .地貌 和 环流 特征 等 进行 综合 判定 。 
(2) 水 文学 上 的 判定 , 即 一 致 区 内 各 站 点 的 极 
值 序列 应 该 为 同一 分 布 ,也 就 是 各 站 点 的 离 差 系 
BN 、 偏 态 系数 和 峰 度 系数 在 一 定 的 容 恳 度 内 保持 一 
致 。 各 站 点 是 否 属于 一 致 区 ,通常 使 用 异 质 性 检验 
指标 也 来 判断 , 当 瓦 <1 时 ,认为 该 一 致 区 可 以 接受 。 
(3) 不 和 谐 检验 ,经 过 气象 学 与 水 文学 上 的 判 
定 后 ,水 文 气象 一 致 区 可 初步 确定 ,然后 计算 出 一 
致 区 内 各 站 点 的 不 和 谐 度 指标 应 ,如 果 不 和 谐 度 旋 
均 小 于 临界 值 , 则 认为 该 一 致 区 划 和 定 合理 ,如 果 存 
在 个 别 站 点 的 不 和 谐 度 D; 超 过 临界 值 , 则 需要 将 不 
和 谐 的 站 点 移 到 周围 区 域 ,甚至 将 该 站 剔除 。 


(GEV) .广义 正 态 分 布 (CCGNO) .广义 帆 雷 托 分 布 
CCGPA) 和 皮尔 逊 亚 型 分 布 (PE3) 最 为 常用 ” ,所 以 
本 文选 用 上 述 5 种 分 布 进行 陕西 省 极端 降水 的 模 
拟 。 首 先 利用 蒙 托 卡 洛 (Monte-Carlo) 模 拟 各 一 致 区 
极端 降水 的 最 优 分 布线 型 ,其 次 ,通过 对 比 一 致 区 
内 各 站 点 的 平均 线性 峰 度 系数 和 所 选择 分 布线 型 
峰 度 系数 之 间 的 差异 ,判断 拟 合 效果 好 坏 ,实际 工 
作 中 ,通常 借助 统计 量 Z"™ 来 判断 , 当 IZ”™|< 1.64 
时 ,认为 拟 合 结果 可 以 接受 ;而 IZ”™| 越 趋 近 于 0, 则 
拟 合 结果 越 好 。 这 里 之 所 以 选择 1.64, 是 因为 Z”™ 
服从 标准 正 态 分 布 , 并 且 以 显著 性 水 平 a=0.1 
确定 。 


2 结果 分 析 


2.1 数据 统计 检验 

地 区 线性 和 矩 方法 要 求 对 数据 序列 进行 充分 的 
统计 检验 ,如 平稳 性 和 独立 性 等 ,以 保证 频率 估计 
值 的 可 靠 性 ,实际 上 ,平稳 性 和 独立 性 也 是 水 文 频 
率 分 析 中 的 2 个 基本 假定 。Mann-Kendall 趋势 检 
验 法 可 检验 数据 序列 的 平稳 性 ,而 一 阶 自 相关 系数 
可 检验 数据 序列 的 独立 性 。 本 文采 用 上 述 2 种 方法 
分 别 对 陕西 省 58 个 气象 站 点 的 极端 降水 序列 
(AMP1、AMP3 、AMP5 .AMP7) 进 行 检验 (图 2)。 

由 图 2a 可 见 , 在 陕西 省 58 个 气象 站 点 的 Mann- 
Kendall 趋势 检验 中 ,AMP1、AMP3 、AMP5 、AMP7 序 
列 分 别 有 32 .28 .33 .35 个 站 点 旦 上升 趋势 ,上 升 站 
点 多 位 于 安康 北部 ,咸阳 宝鸡 北部 .榆林 东部 和 西 
南 角 ,其 中 榆林 西南 角 的 定 边 站 上 升 趋势 最 明显 ， 
而 安康 市 旬 阳 站 的 下 降 趋 势 最 明显 。 所 有 站 点 中 ， 
仅 5 个 站 点 的 个 别 极端 降水 序列 上 升 或 下 降 趋 势 通 
过 了 0.05 的 显著 性 水 平 检验 , 占 全 部 序列 的 4% ,而 
其 余 气 象 站 点 的 极端 降水 序列 均 无 显著 趋势 ， 
此 ,可 以 认为 陕西 省 的 极端 降水 序列 是 平稳 的 。 由 
图 2b 可 见 , 对 于 一 阶 自 相关 系数 , 仅 2 个 站 点 的 
AMPI 序列 通过 了 0.05 的 显著 性 水 平 检验 , 占 全 部 
极端 降水 序列 的 1% ,而 其 他 站 点 的 极端 降水 序列 
均 无 显著 的 自 相 关 , 因 此 ,可 以 认为 陕西 省 的 极端 
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注 : 星 号 标记 表明 通过 0.05 显著 性 水 平 检验 。 
图 2 陕西 省 58 个 气象 站 点 极端 降水 的 Mann-Kendall 趋 热 (a) 和 一 阶 自 相 关 (b) 检 验 结 


Fig. 2 Mann-kendall trend (a) and first-order autocorrelation (b) test results of extreme precipitation 


at 58 meteorological stations in Shaanxi Province 


降水 序列 是 独立 的 。 
2.2 水 文 气象 一 致 区 的 划分 

水 文 气象 一 致 区 的 划分 过 程 , 首 先 , 采 用 开 - 均 
值 法 (依据 站 点 经 纬度 ,高程 和 年 降水 量 ) 进 行 聚 类 
分 析 , 将 陕西 省 划分 为 若干 个 区 ;然后 ,考虑 到 陕西 
省 独特 的 气候 与 地 形 特征 ,结合 不 和 谐 检验 指标 疡 
和 异 质 性 检验 指标 本 ,对 聚 类 分 析 所 得 分 区 进行 适 
当 的 调整 ,使 各 区 内 站 点 尽量 位 于 同一 气候 与 地 形 
条 件 下 ,最 后 ,对 调整 后 的 各 一 致 区 极端 降水 序列 
进行 异 质 性 检验 与 不 和 谐 检验 。 在 调整 过 程 中 , 如 
果 有 分 区 未 通过 上 述 2 个 检验 , 则 将 不 和 谐 的 站 点 移 
到 相 邻 分 区 ,或 者 将 该 分 区 划分 为 2 个 分 区 ,甚至 将 
该 站 点 删除 ,如 此 对 研究 区 的 各 分 区 进行 反复 调整 。 

经 过 多 次 调整 ,本 文 最 终 将 陕西 省 划分 为 6 个 
水 文 气象 一 致 区 , 表 1 为 各 一 致 区 内 站 点 的 异 质 性 
检验 到 结果 。 由 表 1 可 知 ,6 个 分 区 的 各 极端 降水 
序列 均 通 过 了 异 质 性 检验 , 即 甩 均 小 于 1; 同 时 , 除 
了 商 南 站 AMP1 序 列 外 (有 为 1.93 ,略微 超过 临界 值 
1.92) ,所 有 站 点 的 极端 降水 序列 均 通 过 不 和 谐 检 
验 。 图 3 为 陕西 省 6 个 水 文 气 象 一 致 区 最 终 的 划分 
结果 ,与 图 1 对 比 可 见 ,一 致 区 划分 与 陕西 省 的 地 形 


表 1 各 一 致 区 内 气象 站 点 的 异 质 性 检验 用, 结果 
Tab. 1 Heterogeneity test (H) results of meteorological 


stations within each homogeneous region 


一 致 区 AMPI AMP3 AMP5 AMP7 
1 区 0.88 0.19 -1.53 -0.97 
2 区 -1.07 -0.58 -1.34 -1.57 
3 区 -0.44 -2.32 -3.03 -2.84 
4 区 0.98 -1.48 -0.47 -0.96 
5K 0.36 -0.39 0.28 0.28 
6 区 0.95 -0.41 0.49 0.11 


特征 基本 一 致 。 图 3 中 ,1 区 位 于 最 北部 ,有 12 个 站 
点 ,地 形 以 高 原 为 主 ;2 区 位 于 黄土 高 原 , 有 10 个 站 
点 ,地 形 以 高 原 为 主 ,1 区 的 海拔 低 于 2 区 。3 区 主 
要 位 于 渭河 流域 ,有 12 个 站 点 ,地 形 主要 是 平原 和 
丘陵 ;4 区 主要 位 于 秦岭 山地 ,有 5 个 站 点 ,地 形 以 山 
地 为 主 ,海拔 最 高 ;5 区 位 于 秦岭 南部 的 东 侧 ,有 7 个 
站 点 ,地 形 主要 是 丘陵 和 山地 ;6 区 位 于 秦岭 南部 的 
西 侧 , 有 12 个 站 点 ,地 形 也 以 丘陵 和 山地 为 主 。 
2.3 最 优 分 布线 型 的 确定 

考虑 到 AMP1、AMP3、AMP5 和 AMP7 各 序列 的 
结果 比较 类 似 , 下 文 主要 以 AMP1 和 AMP7 为 代表 
给 出 结果 。 表 2 为 各 一 致 区 备 选 分 布 的 拟 合 优 度 检 
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注 : 左 上 和 角 区 域 为 各 一 致 区 及 对 应 一 致 区 AMP1 的 最 优 分 布 。 
图 3 陕西 省 基于 极端 降水 的 水 文 气象 一 致 区 划分 结果 


Fig. 3 Results ofhydrometeorological homogeneous 


region classification based on extreme precipitation 


in Shaanxi Province 


验 结果 。 其 中 ,GEYV 分 布 的 拟 合 效果 最 好 ,GLO、 
GNO 和 PE3 的 拟 合 效果 也 较 好 (对 于 4 种 极端 降水 
序列 ,GEV 共 有 12 次 被 选 为 最 优 的 分 布线 型 ,GLO、 
GNO 和 PE3 则 分 别 有 13、22、13 次 被 选 为 可 接受 的 
分 布线 型 ) 。GPA 的 拟 合 效果 较 差 ,但 在 5 区 的 
AMP5 和 AMP7 序 列 中 , 却 是 唯一 的 可 接受 分 布线 
型 ,其 中 的 原因 还 有 待 进一步 分 析 。 

本 研究 中 数据 序列 长 度 为 45 a, 故 极 值 序列 的 
最 大 值 可 近似 看 作 50 a 一 遇 的 真 值 , 将 其 与 频率 估 
计 值 相 比 较 , 可 以 间接 反映 估计 值 的 准确 性 。 在 前 
文 的 分 区 中 ,1 区 .3 区 .6 区 分 别 位 于 陕 北 、 陕 中 和 陕 
南 ,分 别 属 于 黄土 高 原 .关中 平原 和 秦 巴 山地 ,可 以 
反映 陕西 不 同 地 表 形 态 下 的 降水 特征 。 表 3 分 别 给 
出 这 3 个 分 区 的 AMP1 序列 实测 值 和 50 a 一 遇 的 频 
率 估 计 值 ,以 及 两 者 间 的 相对 误差 (RE)。 可 以 发 
现 ,相对 误差 主要 集中 在 0~30% 之 间 , 有 少数 站 点 
的 相对 误差 超过 30%。 其 中 , 宁 陕 站 的 相对 误差 达 
到 52.4%, 为 相对 误差 的 最 大 值 ,其 实测 最 大 值 超过 
多 年 平均 值 的 4 倍 ,而 相对 误差 较 小 的 站 点 ,实测 最 
大 值 一 般 在 多 年 平均 值 的 2 倍 左右 。 
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a 气象 站 点 N 表 2 各 一 致 区 备 选 分 布 的 拟 合 优 度 检验 结果 
人 Tab. 2 The goodness-of-fit test results of alternative 
国 2 区 GEV g 
= As ae distributions in each homogeneous region 
=o 序列 分 区 O BARA Zw TEAN 
AMPI1 1 GEV 0.29 GEV GNO PE3 
2 GEV -0.67 GLO GEV GNO 
3 GLO 0.55 GLO GEV GNO 
4 GEV -0.56 GLO GEV GNO PE3 
5 GNO —0.09 GEV GNO PE3 
6 GEV 0.05 GEV GNO PE3 
AMP7 1 GEV —0.02 GEV GNO PE3 
2 GEV -0.66 GLO GEV GNO 
3 GEV -0.31 GLO GEV GNO 
4 GEV -0.26 GLO GEV GNO PE3 
5 GPA -0.54 GPA 
6 PE3 0.11 GEV GNO PE3 


利用 绘 点 公式 ,计算 各 站 点 AMP1 序 列 实测 值 
对 应 的 经 验 频率 ,再 确定 该 经 验 频率 下 的 地 区 分 析 
法 估计 值 ,将 两 者 进行 比较 ,可 以 间接 反映 估计 值 
的 准确 性 。 采 用 的 绘 点 公式 如 下 : 

P=(i+A)/(n+ B) (6) 

式 中 ;i 为 各 站 点 AMP1 序列 从 小 到 大 排序 的 序号 ; 
n 为 各 站 点 AMP1 序列 的 长 度 ; 4 和 B 为 绘 点 公式 
参数 ,本 文 4 取 值 为 -0.35, B PUEN 0, 

基于 公式 (6) , 求 出 各 站 点 AMP1 序 列 实测 值 所 
对 应 的 经 验 频率 ,并 基于 最 优 分 布 求 出 该 经 验 频 率 
对 应 的 频率 估计 值 。 图 4 为 各 一 致 区 实测 值 和 频率 
估计 值 的 散 点 图 ,可 看 出 实测 值 与 频率 估计 值 存 在 
较 好 的 线性 相关 ,判决 系数 均 超 过 0.76 , 均 达 到 0.01 
的 显著 性 水 平 。 另 外 ,在 实测 值 较 小 时 , 即 重 现 期 
较 短 时 ,频率 估计 值 较 准 确 ( 散 点 位 置 接近 图 中 的 
对 角 线 ) ,其 中 ,1~3 区 的 实测 值 在 50 mm 左右 时 ,4~ 
6 区 的 实测 值 在 75 mm 左右 时 ,所 对 应 的 频率 估计 值 
最 准确 。 当 实测 值 较 大 , 即 重 现 期 较 长 时 ,频率 估计 
值 明显 小 于 实测 值 , 且 误差 随 实测 值 的 增 大 而 增 大 。 
2.4 单 站 分 析 法 和 地 区 分 析 法 的 比较 

在 单 站 分 析 中 ,各 站 点 极端 降水 的 最 优 拟 合 分 
布 可 由 极 值 图 确定 ,以 陕 北 的 子 洲 站 和 陕 南 的 山 阳 
站 为 例 (图 5)。 由 图 5 可 以 看 出 ,5 个 常用 的 分 布 中 ， 
PE3 分 布 与 子 洲 站 的 AMP1 序 列 拟 合 效果 最 好 ,而 
GPA 分 布 与 山 阳 站 的 AMP1 序列 拟 合 效果 最 好 , 即 
子 洲 站 和 山 阳 站 AMP1 序列 的 最 优 拟 合 分 布 分 别 为 
PE3 和 GPA。 同 理 可 得 其 他 各 站 点 的 最 优 拟 合 分 布 。 
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表 3 研究 区 站 点 AMP1 序 列 的 实测 最 大 值 与 频率 估计 
值 (50 a 一 遇 ) 的 比较 
Tab.3 Comparison of the measured maximum 
value of AMP1 in the study area 
with the estimated quantiles (50-year return period) 


分 区 站 点 Obs(mean) Obs(max) Est-50y RE/% 
/mm /mm /mm 

1 区 府 谷 52.7 181.8 117.0 35.6 
榆林 51.9 105.7 115.3 9.1 
神木 56.8 141.1 126.2 10.6 
佳 县 53.4 216.4 118.5 45.2 
横山 45.2 92.3 100.4 8.7 
子 洲 61.3 150 136.0 9.3 
RAE 57.4 112.5 127.4 13.3 
清润 56.2 98.5 124.8 26.7 
延 川 59.1 157.4 131.1 16.7 
延长 56.3 115.1 124.9 8.5 
宜 川 56.4 113.2 125.1 10.5 
靖 边 47.0 113.2 104.3 7.9 

3 区 ME 58.7 214.6 124.1 42.2 
宝鸡 59.1 169.7 124.9 26.4 
千 阳 63.9 175.1 135.0 22.9 
REUE 53.5 152.4 113.1 25.8 
il 59.3 116.1 125.3 7.9 
白水 65.3 174.2 137.9 20.8 
合 阳 57.5 114.7 121.6 6.0 
乾 县 53.0 117.3 111.9 4.6 
HH 52.6 96.5 111.1 15.1 
兴 平 52.5 146.7 111.0 24.3 
渭南 48.6 102.8 102.7 0.1 
秦 都 46.5 158.5 98.3 38.0 

6 区 略 阳 73.9 161.8 143.8 11.2 
pees 68.5 147.5 133.3 9.6 
城 固 63.1 103.8 122.8 18.3 
TIR 88.2 184.6 171.7 7.0 
南 郑 73.2 151.8 142.5 6.1 
TPE 74.5 304.5 145.0 52.4 
紫阳 93.5 210.8 182.1 13.6 
石泉 77.4 126 150.7 19.6 
DUA 80.8 163.7 157.3 3.9 
镇 巴 121.9 238.2 237.3 0.4 
平 利 69.6 112.4 135.4 20.5 
镇 坪 75.6 174.2 147.2 15.5 


i : Obs(mean) fz AMP1 序列 的 平均 值 ,Obs(max)) 是 AMP1 序列 的 最 


大 值 ,Est-50y 是 50 a 一 遇 的 频率 估计 值 。 


采用 单 站 分 析 法 和 地 区 分 析 法 分 别 进行 计算 ， 
得 到 陕西 省 58 个 站 点 在 不 同 重 现 期 下 2 种 方法 的 
频率 估计 值 ,并 分 别 计算 相对 均 方 根 误差 (RMSE ) 
值 。 相 对 RMSE 越 小 ,说 明 拟 合 效果 越 好 ,准确 性 越 
高 ,图 6 为 AMP1 和 AMP7 序 列 在 50 a 一 过 之 下 相对 
RMSE 值 的 箱 线 图 。 可 以 发 现 ,地 区 分 析 法 的 相对 
RMSE 值 波动 范围 显著 小 于 单 站 分 析 法 的 相对 
RMSE 值 (地 区 分 析 法 的 25% 和 75% 分 位 数 十 分 接 
近 , 因 此 , 箱 线 图 中 的 盒子 很 小 ) ,说 明 地 区 分 析 法 
的 相对 RMSE 值 很 稳定 ,此 结果 与 刘 梦 洋 等 "的 结 
果 基 本 一 致 。 

从 箱 线 图 来 看 ,对 于 AMP1 序 列 , 地 区 分 析 法 相 
XT RMSE 值 大 多 位 于 单 站 分 析 法 相对 RMSE 值 的 波 
动 范围 内 ,两 者 平均 值 的 差距 较 小 ( 仅 2 区 例外 ); 对 
于 AMP7 序 列 ,地 区 分 析 法 的 相对 RMSE 值 较 单 站 
分 析 法 明显 偏 小 ,两 者 平均 值 差 距 较 大 ( 仅 6 区 例 
外 )。 因 此 ,地 区 分 析 法 在 计算 较 长 时 段 极端 降水 
的 频率 估计 值 时 ,准确 性 更 高 。 
2.5 地 区 增长 因子 的 分 布 

一 致 区 内 各 站 点 共有 的 极端 降水 特性 的 频率 
因子 , 称 为 地 区 增长 因子 ,其 实质 是 具有 相同 水 文 
气象 条 件 的 一 致 区 内 ,一定 重 现 期 的 极端 降水 相对 
于 多 年 平均 极端 降水 的 增幅 大 小 。 表 4 给 出 各 一 致 
区 不 同 重 现 期 下 AMP1 序列 的 地 区 增长 因子 。 可 以 
发 现 , 当 重 现 期 为 2 a 一 遇 时 ,地 区 增长 因子 小 于 1， 
且 4、5.6 区 的 地 区 增长 因子 大 于 1、2、3 区; 当 重 现 
期 达到 5 a 一 遇 时 ,地 区 增长 因子 大 于 1, 且 1.2.3 区 
的 地 区 增长 因子 大 于 4.5.6 区 。 另 外 , 随 着 重 现 期 
的 增长 ,地 区 增长 因子 在 增 大 , 且 陕 南 和 陕 北 的 地 
区 增长 因子 差异 也 在 增 大 ,这 是 造成 不 同一 致 区 之 
间 极 端 降水 频率 估计 值 差异 的 一 个 主要 因素 。 根 
据 公式 (5) , 当 重 现 期 很 小 ,如 2 a 一 遇 时 ,其 频率 佑 
计 值 小 于 多 年 平均 极端 降水 , 而 陕 南极 端 降水 事件 
较 多 ,其 2 a 一 遇 频 率 估 计 值 和 多 年 平均 极端 降水 的 
差异 小 于 陕 北 ,导致 陕 南 的 地 区 增强 因子 大 于 陕 
北 。 当 重 现 期 较 大 时 ,如 大 于 5 a 一 遇 , 频 率 估计 值 
可 能 大 于 多 年 平均 极端 降水 , 且 重 现 期 越 大 ,频率 
估计 值 越 大 ,相应 的 地 区 增长 因子 也 越 大 。 同 时 ， 
陕 北 地 区 一 般 降 水 偏 少 ,因而 多 年 平均 极端 降水 也 
较 小 ,但 某 些 年 份 的 极端 降水 很 高 ,从 而 造成 重 现 
期 较 长 时 陕 北 地 区 的 地 区 增长 因子 大 于 陕 南 地 区 。 

一 致 区 内 各 站 点 的 地 区 增长 因子 是 相同 的 , 即 
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图 4 一 致 区 内 各 站 点 实测 值 与 频率 估计 值 的 结果 比较 


Fig.4 Comparison of measured values and estimated quantiles of different stations in the homogeneous region 
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注 : 横 坐 标的 简化 变量 ,为 不 及 制 概率 (通过 -log[-log( 了 | 转换 得 到 的 简化 变量 形式 。 


图 5 子 洲 站 (a) 和 山 阳 站 (b)AMP1 的 极 值 
Fig. 5 The extreme value plot of AMP1 at Zizhou Station (a) and Shanyang Station (b) 


一 致 区 内 的 极端 降水 差异 是 由 各 站 点 本 地 分 量 的 
差异 造成 。 例 如 ,50 a 一 遇 的 情况 下 ,位 于 1 区 的 子 
洲 站 ,其 AMP1 本 地 分 量 达 61.27 mm ,而 1 区 其 他 站 
点 的 AMP1 本 地 分 量 大 多 为 45~55 mm ,小 于 子 洲 站 
的 极端 降水 本 地 分 量 , 这 就 导致 子 洲 站 50 a 一 遇 的 


AMP1 达到 136.01 mm ,而 同 频率 下 1 区 其 他 站 点 的 
AMP1 为 110~125 mm。 另 外 ,不 同一 致 区 的 地 区 增 
长 因子 是 不 同 的 , 即 具 有 相同 的 本 地 分 量 的 情况 
下 ,如 果 站 点 位 于 不 同 的 一 致 区 ,其 极端 降水 频率 
估计 值 是 有 差异 的 。 例 如 ,50 a 一 遇 的 情况 下 , 丹 凤 
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图 6 单 站 分 析 法 和 地 区 分 析 法 频率 估计 值 (50 a 一 遇 ) 的 相对 RMSE 箱 线 图 


Fig.6 Boxplot ofrelative RMSE of the quantiles for single station analysis and regional analysis (50-year return period) 


表 4 一 致 区 不 同 重 现 期 下 AMP1 的 地 区 增长 因子 


Tab. 4 Regional growth factors of AMP1 at different return periods in the homogeneous region 


分 区 2a Sa 10a 20a 50a 100a 200a 500a 1000a 
1 0.91 1.28 1.55 1.83 2.22, 2.54 2.89 3.38 3.79 
2 0.91 1.25 1.49 1.74 2.10 2.39 2.70 315 3.53 
3 0.91 1.25 1.49 1.75 2.11 2.41 2.74 3.20 3.59 
4 0.93 1:25 1.47 1.69 2.00 2.24 2.48 2.83 3.10 
5 0.93 1.23 1.44 1.65 1.93 2.16 2.39 2.72 2.98 
6 0.93 1.24 1.45 1.66 1.95 217 2.40 2.71 2.96 


站 的 AMP1 本 地 分 量 为 61.28 mm, 与 子 洲 站 AMP1 的 空间 分 析 功 能 绘制 等 值 线 图 (图 7)。 由 图 7 可 知 ， 


本 地 分 量 接 近 , 但 由 于 商 南 站 位 于 5 区 ,地 区 增长 因 
子 为 1.93, 而 子 洲 站 位 于 1 区 ,地 区 增长 因子 为 
2.22, 最终 丹 凤 站 的 AMP1 频率 估计 值 为 118.46 
mm ,小 于 子 洲 站 的 136.01 mmo 
2.6 频率 估计 值 的 时 空 分 布 

利用 地 区 线性 矩 法 计算 出 陕西 省 50 a 一 遇 和 
100 a 一 遇 极 端 降水 的 频率 估计 值 ,并 借助 ArcGIS 


4 个 序列 (AMP1 .AMP3 .AMP5 .AMP7) 在 2 个 重 现 期 
下 的 极端 降水 频率 估计 值 空 间 分 布 较为 相似 , 均 为 
南部 极端 降水 较 大 ,而 中 部 的 咸阳 至 商洛 地 区 \ 西 
部 的 延安 西北 角 和 榆林 西部 降水 较 小 ,陕西 省 东部 
的 极端 降水 介 于 前 两 者 之 间 ,总 体 上 大 致 是 从 南 往 
北 减 少 。 基 中 ,6 区 的 频率 估计 值 最 大 ,最 大 值 位 于 
6 区 的 南部 中 心 , AMPI 可 达 270 mm, AMP7 可 达 


频率 估计 值 mm (a) 


频率 估计 值 /mm (c) 频率 估计 值 /mm (d) 
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人 


>z 


图 7 陕西 省 50 a—i8 AMPI (a) ll AMP7(b ) All 100 a 一 遇 AMP1(c) 和 AMP7(d) 极 端 降水 频率 估计 值 的 空间 分 布 
Fig.7 Spatial distribution of the estimated quantiles at 50-year AMP1(a) and AMP7(b) and 100-year AMP1(c) and 


AMP7(d) return period of extreme precipitation in Shaanxi Province 


202109.00044v1 


chinaXiv 


500 mm。 而 在 西 侧 的 黄土 高 原 区 ,AMP1 最 低 只 有 
90 mm, AMP7 U| RÆ 140 mm。 可 见 , 陕 西 省 的 极端 
降水 频率 值 分 布 非常 不 均匀 ,这 与 陕西 省 独特 的 地 
理 分 布 有 关 ( 图 1)。 陕 西 省 南部 海拔 较 低 ,可 接受 
到 南方 输送 来 的 水 汽 ,而 由 于 秦岭 海拔 较 高 (陕西 
省 内 约 2000 m) , 横 跨 陕西 省 中 南部 ,阻挡 了 南方 水 
汽 的 北 进 29 ,导致 中 部 和 和 北部 的 极端 降水 频率 估 
计 值 明显 小 于 南部 。 位 于 秦岭 北 侧 背 风 坡 的 陕西 
省 中 部 地 区 海拔 较 高 的 西部 黄土 高 原 地 区 ,极端 
降水 的 频率 估计 值 最 小 。 将 本 文 结果 与 王 雯 燕 "" 
的 研究 进行 对 比 ,可 以 发 现 ,本 文中 陕西 省 50 a 一 过 
和 100 a 一 遇 极 端 降水 的 频率 佑 计 值 分 布 图 与 王 走 
燕 "" 的 陕西 省 近 60 a 来 极端 降水 强度 、 阅 值 的 分 布 
特征 类 似 , 仅 在 3 区 部 分 站 点 出 现 了 频率 估计 值 偏 
小 的 现象 ,这 可 能 是 由 于 3 区 降水 偏 少 ,而 地 区 线性 
和 抢 法 剔除 了 极端 降水 个 别 异 常 值 而 造成 的 。 总 体 
来 看 ,基于 地 区 线性 矩 法 计算 的 陕西 省 极端 降水 频 
率 估计 值 是 可 靠 的 。 

陕西 省 地 形 多 为 山地 .丘陵 和 高 原 ,极端 降水 
的 发 生 很 容易 诱发 泥石流 和 山体 滑坡 等 灾害 ,即使 
在 极端 降水 频率 估计 值 较 小 的 中 部 和 北部 ,由 于 其 
地 质 条 件 脆 弱 , 突 发 的 极端 降水 也 能 造成 很 大 的 危 
害 。 例 如 ,2010 年 商洛 市 山 阳 县 (属于 5 区) 发 生 的 
“7.23” 特 大 暴雨 ,局 部 地 区 的 24 上 降水 量 超过 190 
mm ,造成 11680 栋 房屋 倒塌 ,15 人 死亡 ;2017 年 7 月 
26 日 ,榆林 市 北部 发 生 特大 暴雨 ,其 中 横山 、 米 脂 、 
子 洲 ( 属 于 1 区 ) 日 降水 量 达 111.1 mm, 140.3 mm, 
218.7 mm, 均 突破 历史 极 值 ,而 强 降水 又 引起 洪涝 、 
泥石流 等 危及 生命 和 财产 安全 的 重大 灾害 。 地 区 
线性 矩 法 具有 对 特大 值 的 稳健 性 ,个别 极 端的 异常 
值 不 会 影响 地 区 频率 曲线 ,计算 单 站 频率 估计 值 的 
时 候 , 是 采用 历年 最 大 值 的 平均 与 地 区 增长 因子 相 
乘 而 得 到 的 。 所 以 陕 北 地 区 出 现 的 极端 降水 异常 
值 ,对 图 7 的 结果 影响 不 大 。 综 上 分 析 ,对 于 陕西 省 
极端 降水 造成 的 多 种 水 文 气象 灾害 ,人 们 应 该 给 与 
更 多 的 关注 。 


3 结论 与 讨论 


本 文采 用 1971 一 2015 年 陕西 省 无 缺 测 情况 分 
布 较 均 匀 的 58 个 气象 站 点 的 逐日 降水 数据 ,运用 地 
区 线性 矩 法 分 析 该 省 极端 降水 的 时 间 和 空间 特征 ， 
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取得 主要 结论 如 下 : 

(1) 研究 区 绝 大 多 数 站 点 的 极端 降水 无 显著 趋 
热 , 也 无 显著 的 自 相关 ,因此 数据 序列 通过 水 文 频 
率 分 析 的 平稳 性 和 独立 性 2 个 主要 的 基本 假设 。 

(2) 利用 K- 均 值 法 对 气象 站 点 进行 聚 类 分 析 ， 
并 结合 不 和 谐 检 验 指 标 D; 和 异 质 性 检验 指标 五, 对 
分 区 进行 适当 调整 ,最 终 将 陕西 省 划分 为 6 个 水 文 
气象 一 致 区 ,其 中 ,GEV 分 布 的 模拟 效果 最 好 ,各 区 
最 优 频 率 估计 值 与 同 频率 的 实测 值 较为 吻合 。 

(3) 相对 均 方 根 误差 RMSE 证 实地 区 分 析 法 的 
优势 ,在 计算 极端 降水 的 频率 估计 值 时 ,地 区 分 析 
法 所 得 结果 较 单 站 分 析 法 更 为 稳健 , 且 对 于 较 长 时 
段 的 频率 估计 值 , 地 区 分 析 法 的 准确 性 更 高 。 

(4) 地 区 分 析 法 的 结果 , 当 重 现 期 为 2 a 一 遇 
时 , 陕 南 的 地 区 增长 因子 大 于 陕 北 ; 当 重 现 期 达到 
5 a 一 遇 时 则 相反 , 且 随 着 重 现 期 的 增长 ,地 区 增长 
因子 在 增 大 ,在 陕 南 、 陕 北 的 差异 也 在 增 大 ,这 也 是 
造成 不 同一 致 区 之 间 极 端 降水 频率 估计 值 差异 的 
一 个 主要 因素 。 

(5) 在 100 a 和 50 a 一 遇 的 重 现 期 下 ,陕西 省 南 
部 的 极端 降水 较 大 ,中 部 的 咸阳 至 商洛 地 区 .西部 
的 延安 西北 角 和 榆林 西部 较 小 ,东部 居中 。 最 大 值 
位 于 6 区 的 南部 ,AMP1 可 达 270 mm, AMP7 可 达 
500 mm , 而 中 部 和 西 侧 的 黄土 高 原 极 端 降水 最 小 ， 
AMP1 只 有 90 mm,AMP7 则 只 有 140 mm。 陕 西 省 极 
端 降水 的 分 布 特征 与 其 独特 的 地 理 特征 有 关 ,尤其 
是 秦岭 的 东西 分 布 阻 挡 了 水 汽 向 北 输送 。 

陕西 省 从 南 至 北纬 度 跨 度 大 ,地 形 地 貌 和 气候 
等 差异 明显 ,极端 降水 在 陕西 省 不 同 区 域 有 较 大 差 
异 。 基 于 地 区 线性 矩 法 对 陕西 省 进行 水 文 气 象 一 
致 区 的 划分 ,这 非常 契合 陕西 省 地 理 特征 空间 差异 
大 的 特征 。 同 时 ,陕西 省 南部 暴雨 多 发 ,北部 黄土 
高 原 区 地 质地 貌 条 件 脆弱 ,极端 降水 的 频率 分 析 对 
其 极端 降水 和 滑坡 、 泥 石 流 等 灾害 的 防 灾 减 灾 有 具有 
重要 意义 。 但 本 次 研究 也 发 现 , AMP1、AMP3、 
AMPS 和 AMP7 各 序列 的 极端 降水 频率 估计 值 在 陕 
西 省 的 分 布 特征 相似 。 这 是 因为 陕西 省 极端 降水 
具有 季节 性 ,多 发 生 在 夏季 ,因此 各 序列 极端 降水 
在 发 生 时 间 上 有 重合 。 因 此 ,下 一 步 作者 将 针对 
AMP1 、.AMP3 .AMP5 和 AMP7 各 序列 极端 降水 频率 
估计 值 的 差异 进行 深入 研究 ,探讨 不 同时 段 对 陕西 
省 极端 降水 的 影响 。 
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Abstract: Extreme precipitation can cause severe disasters in arid and semi-arid regions, such as in the Shaanxi 
Province of Northern China. To investigate the spatiotemporal characteristics of extreme precipitation in the 
Shaanxi Province, this study adopted the daily precipitation data of 58 meteorological stations from 1971 to 2015 
with no missing observations and relatively uniform distribution and used the maximum precipitation on the Ist, 
3rd, Sth, and 7th day to represent the extreme precipitation. The regional L-moments method was further applied 
to study the temporal and spatial characteristics of the extreme regional precipitation, which involves the 
screening and processing of precipitation data, identification of homogeneous regions, goodness-of-fit test, quantile 
estimation for each region, comparison between at-site and regional estimation, deriving regional growth factors, 
and mapping of the spatial patterns of extreme precipitation. The results and conclusions of the study were: (1) 
Shaanxi Province can be divided into six hydrometeorological homogeneous regions, among which GEV 
distribution in each homogeneous region has the best simulation effect, and the estimated optimal quantiles of 
each homogeneous region are in good agreement with the measured value of the same frequency. (2) The 
estimated extreme precipitation quantiles calculated by the regional analysis method have better robustness and 
accuracy compared with the single- station analysis method, especially more significant in calculating extreme 
precipitation over a long period. (3) When the return period is once every 2 years, the regional growth factor of 
southern Shaanxi is greater than that of northern Shaanxi; when the return period is once every 5 years, the 
opposite is true, and with the increase of the return period, the regional growth factor and the difference between 
southern Shaanxi and northern Shaanxi also increase. (4) In the 100-and 50-year return periods, extreme 
precipitation is large in the south, centered in the east, the Xianyang-Shangluo region in the middle, the northwest 
corner of Yan’ an in the west, and the west of Yulin are small. The distribution characteristics of extreme 
precipitation are related to the unique geographical characteristics of Shaanxi Province, especially the east-west 
Qinling Mountains, which block the water vapor transmission to the north, causing differences in extreme 
precipitation between the north and south. 
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